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Abstract 

Power integrity is becoming increasingly relevant due to increases in device functionality and 
switching speeds along with reduced operating voltage. Large current spikes at the device terminals 
result in electromagnetic disturbances which can establish resonant patterns affecting the operation of 
the whole system. 

These effects have been examined using a finite difference time domain approach to solve Maxwell's 
equations for the PCB power and ground plane configuration. The simulation domain is terminated 
with a uniaxial perfectly matched layer to prevent unwanted reflections. This approach calculates the 
field values as a function of position and time and allows the evolution of the field to be visualized. 

The propagation of a pulse over the ground plane was observed demonstrating the establishment of a 
complex interference pattern between source and reflected wave fronts and then between multiply 
reflected wave fronts. This interference which affects the whole ground plane area could adversely 
affect the operation of any device on the board. These resonant waves persist for a significant time 
after the initial pulse. Examining the FFT of the ground plane electric field response showed 
numerous resonant peaks at frequencies consistent with the expected values assuming the PCB can be 
modelled as a resonant cavity with two electric and four magnetic field boundaries. 



Introduction 

Power and signal integrity are of increasing technological relevance due to continuing improvements 
in speed and functionality of electronic systems and the consequent decreasing operating voltages and 
increasing magnitude of switching current for typical digital and mixed signal devices. The rapid 
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switching of large numbers of transistors in an integrated circuit causes a large current spike at the 
power supply terminals resulting in a local ground/power bounce due to parasitic inductances, along 
with an electromagnetic disturbance to the power planes which may establish a resonant pattern as it 
is reflected from the edges of the board. This electromagnetic disturbance may compromise the power 
integrity of other devices or couple to signal traces and affect signal integrity. Here we investigate the 
effects of board resonance in typical multilayer printed circuit boards (PCBs) by solving Maxwell's 
equations for the electromagnetic field using the finite difference time domain (FDTD) method. The 
(FDTD) method gives the electric and magnetic fields as a function of position and time which can 
then be visualised using animation techniques. 

In contrast to frequency domain solutions, a time domain solution captures the response over the 
bandwidth of the stimulus waveform in a single run and can be more efficient for the large 
bandwidths resulting from typical clock speeds in digital and mixed signal systems. The evolution of 
the field distribution can be visualized as the field values are calculated as a function of position and 
time. Moreover all of the frequency domain information can be extracted via a Fourier transform. 

One disadvantage of the FDTD approach is that it is memory intensive as matrices for each 
dimensional component of the fundamental fields are required for each node of the grid. However 
continuing advances in computer technologies allow the technique to be applied to realistic models 
using a standard personal computer (PC). 

Model 

Figure 1 shows the modelled configuration where the PCB power and ground planes are separated by 
a dielectric material and the PCB is surrounded by an air layer. 




Figure 1 Modelled PCB configuration surrounded by air gap. 



The model is discretised into a Yee grid as shown in figure 2 where the E and H field components are 
spatially offset so that each H field component is surrounded by four E field components and vice 
versa. In this finite difference time domain (FDTD) scheme [1] the E and H field components are also 
offset by half of a time step. The vector field components at each node of the Yee grid are calculated 
using explicit or semi-explicit finite centred difference (central-difference) expressions to calculate 
the space and time derivatives in Maxwell's curl equations based on the previous field vector value, 
the values of the field vectors at neighbouring nodes (and a known source). The computation begins 
with all field components set to zero. 
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Figure 2 Yee grid (Taken from [2]). 



The computational domain is terminated by a uniaxial perfectly matched absorbing boundary layer 
(UPML) [3] whereby plane waves of arbitrary incidence, frequency and polarization are matched at 
the boundary to a lossy material for inhomogeneous, dispersive, anisotropic or even non-linear 
domains. Plane waves propagating from the model space decay exponentially within the UPML. This 
is achieved by defining an artificial anisotropic absorbing material with s and fi tensors as shown in 
equations 1 and 2 (where the subscripts 1 and 2 represent the model space and the UPML 
respectively): 
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This variant of the standard FDTD "leapfrog" algorithm allows a unified treatment of both the model 
space and the UPML boundaries with o w and k w are set to and 1 respectively in the model space. 
Starting from the time -harmonic form of Maxwell's curl equation for H i.e 



dH z 


dHy 




SyS Z 





" 


dy 


dz 




s x 


dH x 


_ dHz 


= jd) 





s x s z 





dz 


dx 




Sy 




dH y 


dH x 










s x s y 


- dx 


dy - 




s z - 



(3) 



which can be re -written in terms of D with D x = £ £ r — E x D y = £ £ r — E Vi D z = £ £ r — E. 



Substituting for s in terms of k and a and applying the inverse Fourier transform yields a system of 
time domain differential equations: 
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Equation 4 can be discretized on the Yee grid giving explicit time -stepping expressions for the 
components of D in the UPML from which the E field components can be calculated in a two step 
process. A similar two step process can be used for the H and B field components. The main 
disadvantage of this approach compared to the standard method is the need to store the D and H 
vectors in addition to E and H which effectively doubles the memory requirements. Although the 
UPML is theoretically reflectionless, numerical artefacts can occur due to the spatial discretization. 
To overcome this, k and a are graded over several cells with the outer boundary a perfect electrical 
conductor. A polynomial grading may be used [3] - for example in the x axis: 

<y x {x) = { x / d ) m a xmax , k x {x) = 1 + (K xmax - l)( x /rf) m ( 5 ) 

Thus a x varies from a value of at x = (the PML boundary) to a value of a xmax at the PEC. 

A MATLAB code developed by Willis and Hagness [4] to model a 3D dielectric region with a UPML 
boundary was modified to model the configuration shown in figure 1 . A key modification was to 
allow an anisotropic grid size to capture the smaller geometry in the z axis without the grid becoming 
excessively large. Note a minimum cell size of one tenth of the shortest wavelength of interest is 
required to avoid dispersion effects [5] i.e. 

Ax,Ay,Az<^ (6) 



In addition the Courant-Freidrichs-Lewy (CFL) stability criterion bounds the time step relative to the 
grid discretization to ensure numerical stability [5] and is given by: 

t < j 1 (7) 

Ax^ Ay^- Az^- 

i 



where the phase velocity v is given by v = 



E n £ r 



The power and ground planes are modelled as perfect electrical conductors (i.e. the tangential E field 
components E x and E y are set to zero throughout). 

Results 

Figure 3 shows the evolution of the z component of the electric field at the ground plane as a function 
of time where the voltage source is located at the centre of the board. In this case the edge of the 
board is terminated by a UPML boundary. 

Excitation was provided by a current source implemented by introducing modified semi -explicit field 
update coefficients at the source location [6]. A differentiated Gaussian pulse excitation is used given 
by equation 8 with J = 1 A/mm 2 , t = 100 pS and x = 50 pS. 

/ t =/o(t-t )e"^) (8) 
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Examination of the Fourier transform of the current pulse confirms it has significant components to 
frequencies > 15 GHz. 

The board dimensions are L x = 100 mm, L y = 60 mm and L z = 0.8 mm and the size of the Yee grid is 
Ax = Ay =2 mm and Az = 0.4 mm. Hence according to equation 6, dispersion should not be evident for 
frequencies up to 30 GHz. A time step of 1 pS was used which is well below the CFL stability limit of 
1.33 pS given by equation 7. 

Adequate convergence using only 2 cells in the z-axis was confirmed by running a smaller model in x 
and y dimensions with 1 to 15 cells in the z-axis where no significant difference the electric field ( < 
0. 1 %) was seen. 



(a) 

Ground Plane E z , time = 2Dps 




(b) 



Ground Plane E z , time = lOOps 
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(e) (f) 



Ground Plane E z , time = 300ps Ground Plane E z , time - 5rj0ps 




Figure 3 Evolution of ground plane E z distribution with UPML boundary. 

From figure 3 it can be seen that the wave front spreads radially over the surface of the board and has 
decayed to zero by ~ 600 ps. Figure 4 shows a similar model with an air layer included between the 
edge of the PCB and the UPML boundary. In this case the wave front is reflected from the 
dielectric/air interface in they-axis after ~ 300 ps. A complex interference pattern is set up between 
source and reflected wave fronts and then between multiply reflected wave fronts which affects the 
whole power and ground plane area and could adversely affect the operation of any device on the 
board. These resonant waves persist for a significant time after the initial pulse. Although it should be 
noted that this is overestimated as the model uses a lossless dielectric and perfect (i.e. lossless) 
electric conductors, although in practice these losses will be small at the frequencies of interest. The 
model does capture the small transmission loss which does occur as there is some field coupling 
across the dielectric/air boundary. This too is of interest as the radiated fields can result in 
electromagnetic interference [7] 
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Figure 4 Evolution of ground plane E. distribution for PCB surrounded by an air layer. 

The power planes can be modelled as a resonant cavity with two perfect electric conductor 
boundaries (corresponding to the power and ground planes) and four perfect magnetic conductor 
boundaries corresponding to the dielectric/air interfaces when the resonant frequencies are given by 
equation 9 [8, 9] 
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where c is the velocity of light, s r is the relative permittivity, m and n are integers representing the 
order of the resonant mode and L x ,L y and L z are the board dimensions. Note that the small addition of 
(LJ2) to the L x , and L y dimensions in equation 8 is intended to approximate the effects of fringing 
fields at the dielectric/air interfaces. 



A fast Fourier transform (FFT) of the ground plane E z did not show the expected resonant behaviour 
presumably because the response is dominated by the source waveform. In addition it has been 
reported that not all modes are excited equally, there is some dependence on the location of the source 
and some modes can be suppressed depending on the location of the source and measurement ports 
[9]. To avoid any possible mode suppression the source was moved a corner of the board and the 
ground plane E~ measured at the opposite corner. The current source as replaced by a resistive voltage 
source implemented using the appropriate semi-explicit update coefficients [6] with amplitude V = 
5Fand source resistance Rs = 10Q. This allows any effects due to source match or mismatch to be 
investigated by varying the source resistance. The voltage excitation used the same differentiated 
Gaussian pulse shape as for the current source. The board size was reduced to L x = 40 mm, L y = 30 
mm and L z = 0.8 mm to increase the resonant frequencies to allow them to be resolved in fewer time 
steps. 

Figure 5 shows a FFT of the ground plane E z measured at the opposite corner to the source after the 
model was run for 3 x 10" 8 s (30000 steps). The resonant peaks are generally well resolved. Table 1 
compares the frequencies of the peaks with the values expected from equation 9. The measured and 
calculated frequencies are in excellent agreement. The modes labelled J and O are only just resolved 
and were assigned by closely examining the spectrum. However it should be noted that that the 
predicted frequency separation from their neighbours is small. Unfortunately it was not possible to 
improve the resolution by increasing the number of time steps further without exceeding the available 
memory of the 64 bit / 8GB PC. 




Frequency (GHz) 

Figure 5 Frequency spectrum of ground plane E. for board size L x = 40 mm, L y = 30 mm and L z = 0.8 mm. 
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3.03 


2 





3.71 
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4.46 
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4.93 
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4.79 
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5.27 
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5.13 


3 





5.57 
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5.44 
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6.09 


I 


6.06 
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6.18 
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6.18 
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7.40 
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7.14 


4 





7.43 
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7.30 


3 


2 


7.44 
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7.38 


1 


3 


7.63 
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7.69 


4 


1 


7.82 


O 


7.75 


2 


3 


8.28 


P 


8.16 


4 


2 


8.92 


Q 


8.74 


3 


3 


9.26 


R 


9.05 


5 





9.28 


S 


9.39 


5 


1 


9.60 


T 


9.52 





4 


9.87 


U 


9.70 


1 


4 


10.04 


V 


10.23 



Table 1 Comparison of measured and calculated resonant frequencies 
for board size L x = 40 mm, L y = 30 mm and L. = 0.8 mm. 

Conclusion 

Power integrity is becoming increasingly relevant due to increases in device functionality and 
switching speeds along with reduced operating voltage. Large current spikes at the device terminals 
result in electromagnetic disturbances which can establish resonant patterns affecting the operation of 
the whole system. 

These effects have been examined using a finite difference time domain approach to solve Maxwell's 
equations for the PCB power and ground plane configuration. The simulation domain is terminated 
with a uniaxial perfectly matched layer to prevent unwanted reflections. This approach calculates the 
field values as a function of position and time and allows the evolution of the field to be visualized. 

The propagation of a pulse over the ground plane was observed demonstrating the establishment of a 
complex interference pattern between source and reflected wave fronts and then between multiply 
reflected wave fronts affecting the whole ground plane area. This wave interference effect could 
adversely affect the operation of any device on the board. These resonant waves persist for a 
significant time after the initial pulse. Examining the FFT of the ground plane electric field response 
showed numerous resonant peaks at frequencies in excellent agreement with the expected values 
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assuming the PCB can be modelled as a resonant cavity with two electric and four magnetic field 
boundaries. 



References 

1 . Taflove, A. "A Perspective on the 40-Year History of FDTD Computational 
Electrodynamics", Applied Computational Electromagnetics Society Journal, 22 (1), 2007. 

2. Kung, F. "Modeling of Electromagnetic Wave Propagation in Printed Circuit Board and 
Related Structures", PhD Thesis Submitted to Multimedia University, Malaysia, 2003 

3. Gedney, S. "Perfectly Matched Layer Absorbing Boundary Conditions", 2005a , in: Taflove, 
A. and Hagness, S. Eds. "Computational Electrodynamics: The Finite Difference Time 
Domain Method", 3 rd Edition, Artech House, MA. 

4. Willis, K.J. Hagness, S.C. "3-D FDTD code with UPML absorbing boundary conditions", 
2005. [online: 

http://www.artechhouse.com/Default.aspx7eAppT ype=l&strFrame=Misc/taflove.html1 

5. Taflove, A. & Hagness, S. "Computational Electrodynamics: The Finite Difference Time 
Domain Method", 3 rd Edition, Artech House, MA, 2005. 

6. Piket-May, S.C, Gwarek, W, Wu, T.L., Houshmand, B., Itoh, T., Simpson, J. " High-Speed 
Electronic Circuits with Active and Non-linear Components", 2005 in: Taflove, A. and 
Hagness, S. Eds. "Computational Electrodynamics: The Finite Difference Time Domain 
Method", 3 rd Edition, Artech House, MA. 

7. Kim, J. Pak, J., Park, J., Kim, H. "Noise generation, coupling, isolation, and EM radiation in 
high-speed package and PCB", IEEE International Symposium on Circuits and Systems, 
2005. 

8. Chen, R.L., Chen, J., Hubing, T.H., Shi, W. "Analytical Model for the Rectangular Power- 
Ground Structure Including Radiation Loss" IEEE Trans. EMC, 47 (1), 2005. 

9. Lei, G.T., Techentin, R.W., Gilbert, B.K "High frequency characterization of power/ground 
plane structures" IEEE Trans. Microwave Theory and Techniques 47 (5), 562, 1999. Van den 
Berghe, S., Olyslager, F., De Zutter, D., De Moerloose, J., Temmerman, W. " Study of the 
Ground Bounce Caused by Power Plane Resonances" IEEE Trans on Electromagnetic 
Compatibility, 40 (2), 111, 1998. 



Page 11 of 11 



